Boson Realization of the SU(4) Model of High- Temperature Superconductivity 
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The SU(4) algebraic model of high-temperature superconductivity is studied employing the boson 
mapping techniques. The bosonization of the model enables us to get exact numerical solution 
of the model. The order parameters are discussed. The situation close to the SO(5) dynamical 
symmetry limit is interpreted as a modelling case for the behaviour of D-wave superconducting and 
antiferromagnetic phases in cuprates. 



The high-temperature superconductivity represents a 
challenging task for the theory of many-electron corre- 
lated systems. Data suggest that it is related to the sin- 
glet D-wave pairing (dSC) phase. One of the intriguing 
aspects is then a proximity of the dSC phase to the an- 
tiferromagnetic (AF) phase. Indeed, in cuprates, the AF 
order is observed at half filling whereas the dSC order 
develops as the material is doped by holes. 

The interplay between the dSC and AF phases is 
treated by employing the symmetry principles in the 
SO(5) model proposed by S.-C. Zhang §. Two dSC or- 
der parameters and three AF order parameters are uni- 
fied in a five dimensional vector. The SO(5) rotations 
induce then transitions between two phases. Microscopic 
models exhibiting the SO (5) symmetry have also been 
discussed 

Recently, Guidry and collaborators have embedded the 
SO(5) algebra into a larger SU(4) algebra The mi- 

croscopic realization of the SU(4) generators has been 
given and the three dynamical symmetry limits of the 
SU(4) model have been studied. The SO(5) limit is in- 
terpreted as a critical dynamical symmetry interpolating 
between the dSC and AF phases. 

In the present paper, we discuss the boson realization 
of the SU(4) algebra. The bosonization of the model en- 
ables us to get exact numerical solution of the model. 
The closed forms for the dSC and AF order parameters 
in the dynamical symmetry limits are obtained. The in- 
terplay between two phases is studied also outside these 
limits. In agreement with Ref. |J, the situation close to 
the SO (5) limit is suggested as a modelling case to the 
real behaviour in cuprates. 

The 15 SU(4) generators are microscopically con- 
structed on a lattice with the degeneracy Q as B 
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Here, c t . denotes the creation operator of a fermion 
with momentum k and spin projection i — Q is 

an AF ordering vector so that exp(iQ • x) = ±1 at every 
lattice site, and g(k) = sgn(cos k x — cos k y ) with the 
properties <?(k + Q) = - 5 (k) and |<?(k)| = 1. §0 

There are three dynamical symmetry chains conserving 
total spin S and particle number n (or charge M = \{n — 
tt)) discussed in {|||: 

i) the SO (4) limit comprising a subset of operators of 
staggered magnetization and spin S 1 ^, 

ii) the SO (5) limit with the triplet 7r-wave creation and 
annihilation operators ttI and 7f M , and S^, 

iii) the SU(2) limit with the singlet D- wave creation and 
annihilation operators and D (and S 1 ^). 

One can move between the dynamical symmetry limits 
with a Hamiltonian 

ffl = -y((l + -(l-x)^ ■%) (2) 

with x = 1 for the SU(2) limit, x = -1 for the SO(5) 
limit, and x = for the SO(4) limit. The SU(4) quadratic 
Casimir operator 

Csu(4) = D^D -^■% + Q-Q + S-S + M(M - 4) 

gets for the most symmetric collective SU(4) subspace 
(which is only discussed in the present paper) the value 
iQ(Q + 8). Then, an alternative parameterization of the 
SU(4) Hamiltonian up to uninteresting for the present 
discussion constant, charge, and quadratic spin terms is 
written as 



-G u ((l-p)D^D+pQ-Q) 



(3) 



with p = for the SU(2) limit, p = 0.5 for the SO(5) 
limit, and p = 1 for the SO (4) limit. 
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Even though the Hamiltonians (g) and (||) must gen- 
erally be equivalent, the assumed range of the parame- 
ters Gi,G u > 0, 1 > x > -1, and 1 > p > leads 
to different physical picture. For the SU(2) limit, both 
parameterization agree. For the SO(4) and SO(5) limits 
and in-between, however, the positive values of the inter- 
action strengths G\ and Gn imply different signs of the 
Hamiltonians (Q) and (|J). The spectra of two parame- 
terization are reversed with the ground state of (Q) being 
the highest-lying state of (^|) and vice versa. 

For the dynamical symmetry limits, the SU(4) Hamil- 
tonian is analytically solvable. Outside the limits, the 
numerical diagonalization must be performed. In both 
these tasks, an application of the boson mapping tech- 
niques proves to be useful. The Dyson boson realiza- 
tion of the SU(4) algebra is constructed by mapping the 
bifermion operators (|l|) onto bosonic operators formed 
from the scalar boson operators d and vector boson op- 
erators p M 
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By employing (||), one can easily construct the boson 
image of the SU(4) Hamiltonian. The bosonized task is 
relatively simple to solve numerically in the boson space 
SU(l) rf ®SU(3) p [§. To obtain the spectrum of S = 
states including the ground state, it is sufficient to diag- 
onalize a matrix of dimension n/2. That represents an 
enormous truncation of the original fermion space. 

In the boson treatment, we obtain the exact and full 
solution of the fermion problem. The eigenstates states 
have good charge and spin quantum numbers which re- 
flect the symmetry of the original task. For them, we 
relate the dSC order parameter to the matrix element of 
the D-pair transfer between the ground states 
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whereas the AF order parameter is connected to the re- 
duced matrix element of the Q operator between the S=0 
ground state and the lowest lying S=l state 

0o = (nS = l\\Q\\nS = 0) . 

Even if the and Q order operators do not always 
belong to the generators of the particular dynamical sym- 
metry limit chains, we have succeeded to obtain the ana- 
lytical formulas by inspecting the numerical results. We 
have got 



SU(2) parameterizations (2) and (||) 
al = -(ra + 2)(2ft-n) 
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SO (5) parameterization 
2 _ (2ft - n)(2ft - n + 2^ 
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SO (5) parameterization (3) 
l^_ (n + 2)(2ft-n + 6) 
3 2ft - n + 8 
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SO (4) parameterization 
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(2ft-n + 6)(n + 2) for n = 4k + 2 



Pi = for n = 4k 



SO(4) 



parameterization (3) 

I(n-n)(n + 4) 
1 



Pa = 7 n ( n + 4 ) 
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The results (§) and |]) for the SO (5) and SO (4) limits 
of the parameterization (||) are remarkable. In the SO (4) 
limit, the AF order parameter /3 increases linearly with 
n up to the value ft/2 for n — ft. The dSC order param- 
eter ao — > when n — > and n — > ft, and reaches the 
maximum value 4ft at the quarter filling n = ft/2. 

In the SO(5) limit, the order parameters for n <C ft 
agree with the values of the SC case (H) since then the 
SO(5) ground state of the parameterization (||) has the 
form of the D-pair condensate |g.s.) oc D> 2 ra |0). The 
situation changes for n — > ft when a steep increase ap- 
pears in the AF order parameter whereas the SC order 
parameter goes to zero. One may get more insight into 
this phase transition by introducing the fractional dop- 
ing of holes (5 = 1 — n/ft. Then for ft — > oo, Eqs.(ph are 
rewritten as 
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from which form the difference between the half filled 
case (5 — 0) and the hole doped situation (5 > 0) is 
easily seen. 
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FIG. 1. AF order parameter (bottom panel) and dSC order 
parameter (top panel) as a function of the fractional doping of 
holes 5. Degeneracy of the laticce f2=2000. Curves are shown 
for the values of the parameter p in (j^J) p=l (SO(4) limit), 
p=0.5 (SO(5) limit) and two intermediate cases p=0.52 and 
p=0.75. 

In Fig.Q, the order parameters are shown as obtained 
with Hamiltonian (0) for £1 = 2000. Besides the dynam- 
ical symmetry limits SO(5) (p—0.5) and SO(4) (p=l), 
two calculations are given for the values of the parame- 
ter p=0.52 and 0.75. Outside the dynamical symmetry 
limits, the order parameters are close to the SO (5) values 
for 8 — ► 1 and to the SU(4) values for 8 — > 0. The re- 
sults for the AF order parameter in the case p=0.52 are 
very much similar to those obtained in Ref. || with the 
coherent state approach. 

We thus confirm the findings of Refs. B and that 
the SU(4) algebraic model in a situation close to its 
SO (5) dynamical symmetry limit is a convenient tool to 
model behaviour of high-temperature-superconducting 
cuprates. The transition from the AF phase to the dSC 
phase naturally appears as the system is doped by holes 
starting from the half-filled case. 

Attractive feature of the SU(4) algebraic model is its 
exact solvability. Then despite the fact that some por- 
tion of the real physics might not be present in the model 
(there is no kinetic-energy term in the SU(4) Hamilto- 
nian), it can still be useful in testing approximate many- 
body procedures. For the dynamical symmetry limits we 



have obtained closed expressions from which behaviour 
of the order parameters in the phase transition region is 
seen explicitely. 
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Our 7r T operator differs by a multiplicative factor i from 
Ti-t of Ref. " 



k k, 



k L k+Q ,i -k, 



! ) and quasi-spin 77-pair operators (ry = 
! ) , besides the D-pair and 7r-pair oper- 



Note, that the SU(4) algebra can be embedded into 
the SO(8) algebra containing the Cooper S-pair (S^ = 

E 
E 

ators. 

[8] The spurious (nonphysical) states arise in the boson 
space for n > Q which is outside the focus of our analysis. 
The Dyson mapping (^) is nonunitary leading to nonher- 
mitian boson Hamiltonians. Standard procedures exist 
how to treat the original hermitian fermion problem in 
the nonunitary boson method (see for example A. Klein 
and E.R. Marshalek, Rev. Mod. Phys. 63, 375 (1991)). 

[9] Of course, an alternative methods for the numerical solu- 
tion may employ the algebraic properties of the fermionic 
SU(4) algebra and any of its dynamical symmetry chains. 
[10] In the dynamical symmetry limits both in the parameter- 
izations (^) and (H), the first excited S=l state exhausts 
fully the Q-strength from the ground state. Similarly, the 
D-pair transfer from the ground state goes fully or almost 
fully to the ground state of the adjacent system. An ex- 
ception is the SO(5) limit in the parameterization ^ for 
n = 4k + 2 case. 
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